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Abstract 

We present some simple scenarios for string breaking on the lattice based 
on a crude strong coupling model introduced previously. We review the de¬ 
pendence of the model on lattice spacing and extend it to include degenerate 
dynamical quarks and also meson exchange diagrams. A comparison is made 
between quenched and unquenched calculations. We examine string breaking 
in the presence of a static quark-diquark system, a situation that is specific to 
SU( 3) . 
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1 Introduction 


String breaking has importance in lattice calculations because it reveals directly the 
effect of dynamical quarks. However the energy crossover effect characteristic of string 
breaking has proved difficult to observe i-0. Recently progress has been made and 
string breaking has been observed as a mixing phenomenon directly in SU(2) Higgs 
models in three and four dimensions i 11- The phenomenon is currently under 
investigation in the much more demanding case of QCD [| 


12 . A crude model of the 


process based on strong coupling ideas has been proposed which pictures it as a mixing 
process between a string state and a two-meson state [I 


14 . It shows that because of 


its dependence on a mixing angle the Wilson loop may not reveal the string breaking 
energy crossover unless the mixing region is sufficiently broad and that in any case it 
is not by itself a satisfactory observable for revealing the crossover. This is consistent 
with numerical simulations j3j] which show the need for a suite of operators with 
members appropriate to both the string and two-meson states. 

In this paper we develop the crude strong coupling model to discuss a number of 
simple mixing scenarios and use it to discuss the relationship between quenched and 
unquenched calculations. We take the opportunity to improve the model by including 
energy factors associated with the emission of light quarks from bound states. These 
factors are required in order to achieve results independent of the lattice spacing. They 
were missed in the original presentation ]H| but are essential parameters Jl4|. [15| . 

We elaborate the model on the one hand by including several quark flavours and on 
the other by including diagrams that represent meson exchange. 

Finally we apply our approach to more complex situations that are specific to SU(3) 
and examine the case of string breaking in the presence of a static quark-diquark 
system. 


2 Rules for Simple Planar Model 


The model we use was explained previously |T3|, |1J]. It is based on simple planar 
lattice diagrams with internal dynamical quark loops. The initial formulation of the 
model is in terms of the rules for a “strong coupling expansion” ||T6|, [TTH, modified to 


accommodate the presence of bound states in a simple way. The rules are (slightly 
modified from those of ref |T3| as indicated in ref ||14|| ): 

1. A factor of e~ a for each plaquette, where a is the (dimensionless) string tension. 
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2. A factor of 


2k 


1 + 7-e 


for each Wilson quark line in the direction of the unit vector e, interior to the 
diagram, ffere k is the standard quark hopping parameter. 


3. A factor of 


2k' 


m 


for each Wilson quark line propagating in the direction of the unit vector e along 
a static anti-quark line. This creates the model for the static-light bound state 
meson. 


4. A factor of \/W for each emission of a light quark from a meson state, ffere W is 
a (dimensionless) energy parameter that determines the rate of emission of light 
quarks from the meson. Such a parameter is essential in order that the results of 
the model are independent of the lattice spacing. 

5. A factor of (—1/3) (in the case of £77(3)) for each internal quark loop. 

6. A trace over the spin matrix factors for each internal quark loop. 


In the above the hopping parameter, k, reflects the light quark mass, 2k = e~ mq , and 
the bound state hopping parameter, k', determines the energy of the static meson, 
2k! = e~ EM . (In ref [[L3| Em was the energy of a two-meson state and therefore had 
twice the value assigned here.) 

If we apply the above rules to the diagram in Fig. I| we obtain for the string-string 
transition amplitude the result 


Gss = (e" CTR ) Tl -±=We~ mqR ( 


„-2 E m \ T2 J_ 

> v/6 


We~ 




t 3 


( 1 ) 


3 Lattice Spacing Dependence 

The string breaking scenario considered previously involved a static quark and anti¬ 
quark separated by a (dimensionless) distance R. In the notation of previous papers 
00, we have 

a = e ~ aR , b = e~ 2EM , c = 4= We~ mqR . (2) 

V6 
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The basic string state has an energy aR and the two-meson state an energy 2 Em- As 
shown in ref [|13j the propagating eigenstates have energies E± = — log A± where 



+ b ± \]{a — b ) 2 + 4afec 2 | 


(3) 


To pick out the lattice spacing independent part we can treat all the dimensionless 
parameters as O(a^), where cll is the lattice spacing. We then treat R and any time 
interval, T, involved as 0{aJ j 1 ) and expand to the lowest significant order in cll- We 
then find 


E± = i | aR + 2 E m ± J{aR - 2 E M ) 2 + ^W 2 e~^ R 


(4) 


Having carried out this operation we can treat the above formula as referring to di- 
mensionfull quantities. Note that it is essential to have the energy parameter W as a 
factor in the quantity c in order that this result holds true [1141 ITT]. 


Similar remarks apply to other quantities in the model. The mixing angle is 


tan 6 = 


-{a — b) + \J(a — b) 2 + Aabc' 2 


2 \fabc 

Carrying out the expansion in small quantities we find 

aR - 2 E m + ^J{aR - 2 E M ) 2 + \W 2 e~ 2m * R 


(5) 


tan 9 = 


| We- m o R 


( 6 ) 


The formulae in ref |H| for the mixing range A R and the energy split at maximal 
mixing, A E are correct when the energy parameter W is included. As indicated in ref 


lj, the results are 


AR = 7T 



A E 



(7) 


The relation 


AR = 7T 


A E 


( 8 ) 


a 

is preserved. All these formulae may now be regarded as being explressed in physical 
parameters. The transition energy parameter is set by A E. Eq (|D is then a prediction 
for the mixing range in terms of the maximal energy split. 
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4 Quenched/Unquenched Comparison 


An important reason for measuring string breaking is to identify the effect of dynamical 
fermions directly. It is therefore useful to consider the comparison between relevant 
quantities in the quenched and unquenched situations. A natural amplitude to measure 
is the transition amplitude between the string and the two-meson state 0. The 
unquenched case was considered previously 0. [14|. The result was an amplitude of 
the form 

Gms(T) = \fab sin 0 cos 6 — A^ -1) ) , (9) 


where T is the full interval in the time direction associated with the measurement. The 
significance of this form is that the amplitude is suppressed outside the mixing region 
where the factor sin 9 cos 6 vanishes. If we reexpress the formula in terms of physical 
parameters then we find 


1 / x [\We- m ^ R 

Gms(T) = - 

2 \\J(crR — 2E m ) 2 + |hh 2 e 

The corresponding quenched approximation can be described in our model as a sum 
over graphs of the form shown in Fig. |j. The resulting amplitude expressed in terms 
of the parameters a, b and c is 

GmI(T) = £ b T *ca T ' , ( 11 ) 

Ti=l 

where T = 7\ + T 2 .The result is 

n(Q)frr\ ,b T ~ 1 -a T - 1 

Gms\T) — cab - — - . ( 12 ) 

When translated into physical parameters we find 





g—2 E m T g— crRT\ 


(13) 


The quenched and unquenched results are very similar in the region \aR — 2 Em\ > 
^/| We~ mqR . In the complementary region, the inner part of the mixing region on our 
definition, they exhibit very different behaviour. The unquenched amplitude exhibits 
the continuous crossover while the quenched amplitude shows an abrupt change from 
one exponential behaviour to the other at the crossover point. It would be extremely 
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interesting if this contrast in behaviour could be verified in a real simulation. Although 
our model rather oversimplifies the relationship of the quenched and unquenched cal¬ 
culations it is still possible that a measurement of the above transition amplitude in 
a quenched simulation would yield an estimate of the energy W that could stand as a 
prediction for the unquenched calculation. 


5 Several Flavours 


So far the model has been formulated only for a single flavour of light quark. The 
presence of several flavours of quark is accounted for by summing over the Nf flavour 
insertions of quark loop. This means that the model has Nf + 1 channels, namely the 
string channel and the Nf two-meson channels associated with the different flavours of 
quark. It is sufficient to consider the case Nf = 2 . Following the same analysis as in 
the single flavour case and assuming that the light quarks are degenerate we easily see 
that the summation of relevant graphs is achieved by considering a 3 x 3 correlation 
matrix G(T) that obeys the update equation 



G(T + 1) = AG(T) , 

(14) 

where 

A = 

(a ac ac \ 
be b 0 , 

^ 6c 0 b J 

(15) 


and we can assume that G{T) = A T . It is immediately obvious that the two-meson 
state that is anti-symmetric under the interchange of flavours propagates with the 
factor b T = ^t' 2EmT . It does not mix with the string state. The symmetric state does 
mix with the string. If we choose the symmetric and antisymmetric states as basis 
states then the matrix A assumes the form 


A 


f a y/2 ac 0 \ 

\[%c b 0 , 

v 0 0 b J 


(16) 


The problem therefore remains a two-channel mixing problem of the same form as 
before with the minor change that c —> \/2c . There is no need to carry the explicit 
calculation further since the results are obvious. 

The general case of Nf flavours follows the same lines. Only the completely flavour 
symmetric combination of two-meson states mixes with the string. The mixing being 
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of a strength \jNf times greater than the single flavour case. The remaining Nf — 1 
two-meson channels orthogonal to the symmetric one remain unmixed. 

If the light quarks are not flavour degenerate then to a first approximation the two- 
meson state corresponding to the lightest quark will show the first and strongest mixing 
with the string state since it lies at the lowest energy. Heavier quarks will mix at larger 
distances with decreasing strength. It would be interesting to test this obvious and 
plausible scenario in a full QCD calculation. 

6 Extended Planar Model 

It was implicit in the summation technique used previously that diagrams of the form 
shown in Fig. [3] were not included [|l|, [14|. To include them requires a new rule for 
the vertex associated with the emission of the light quark pair. It is 

7. A factor \JW P for the pair emission vertex, where W p is a dimensionless energy pa¬ 
rameter. Again this factor is necessary in order to eliminate explicit dependence 
of the results on the lattice spacing. 

If, after emission, such a light quark-anti-quark pair were to propagate along a link, 
the corresponding factor according to the rules would be 



Naively the hopping parameter factor (2k) 2 is equal to e~ 2mq . In the context of the 
model it would not be unreasonable to identify the quark pair exchange as a meson 
exchange. If the meson has a mass m then we should replace (2k) 2 with e~ m . This is 
the interpretation we will adopt. 

For completeness we give the list of rules for computing diagrams that include strings, 
two-meson states and light meson exchanges. 

1. A factor a = e~ aR that propagates the string by one time step. 

2. A factor b = e~ 2EM that propagates the two-meson state by one time step. 

3. A factor c = (IF/a/6) (^k) R associated with the transition from string to two- 
meson state and vice-versa. 

4. A factor 5 = | W p e~ mR for each quark-anti-quark exchange. 


7 




As in the original model the summation can be carried out by computing a 2 x 2 matrix 
of propagators connecting the string and two-meson states 


G{T) = 


Gss(T) Gsm(T ) 

G M s(T) G M m{T) 


which obeys the update equation 


G(T + 1) = AG(T) , 


where 


A = 


a ac 
be 6(1 + <5) 


(18) 


(19) 


( 20 ) 


The natural solution to eq ([19D is G(T) = A T . As before A can be expressed in the 
form A = DOAO~ 1 D~ 1 , where D and O have the forms 


D = 


y/a 0 

0 Vb 


O = 


cos 9 — sin 9 
sin 6 cos 9 


( 21 ) 


and A is a diagonal matrix of eigenvalues, 

1 


A± = - |(a + 6(1 + 6)) ± \J(a — 6(1 + 6)) 2 + 4a6c 2 | . 


( 22 ) 


The mixing angle 9 , that describes the overlap of the string and two-meson channels 
with the eigenchannels of def ini te energy, E± = — log A± , satisfies 


tan# = 


-(a - 6(1 + 6)) + \J (a - 6(1 + 6)) 2 + 4a6c 2 
2\fab c 


(23) 


The analysis of mixing is essentially the same as before except that the centre of the 
mixing region where 9 = 7t/4 occurs when a = 6(1 + 6) that is when R = R c where 


aR c = 2 E m + log(l + 6) ~ 2 E M + -W p e~ mRc 

o 


(24) 


This is consistent with the idea that the extra light meson exchange represented by 
S corresponds to a potential interaction between the two static mesons of the form 
(Wp/6)e~ mR that displaces the centre of the mixing region in the appropriate way. If 


we estimate the mixing region as A R = fyy|/?=/i c then we find 

[2 We~ mqRc 

AR = 7 t\ 


3a — AmW p e~ mRc 


(25) 







This is essentially the same result as derived previously | [TT| . |T4|| but now including the 
effect of the potential between the mesons. As before the energy split at the point of 
maximum mixing is A E = 2c, that is 


A E = 



We have therefore the relation between A R and A E of the form 


(26) 


A R = 7T- 

a — 


A E 

| mW p e~ mRc 


(27) 


These results may now be interpreted as referring to quantities in physical units. In 
deriving them we have dropped corrections that vanish with the lattice spacing. 


7 Quark/Diquark String Breaking 


So far string formation and breaking in QCD has been sought in a context in which 
the string joins a static quark and anti-quark pair. An interesting alternative scenario 
involves a static-light diquark supporting a string along with another static quark. 
This is a situation that is particular to SU(3 ) gauge theory. We again have quenched 
and unquenched possibilities. 

A relevant correlation function can be formulated in the following way. 


T(T) = (Ic/(C 1 )„ lft C/(C2)„ftG,„, ft (T) 

taia2&3 


(28) 


Here the angle brackets indicate averaging over the gauge fields and G qcc3/ 3 3 (T ) is the 
light quark propagator from the origin to the origin over a time interval T . The 
matrices U(Ci ) and U(Ci) are each the product of gauge SU(3 ) matrices along the 
paths C\ and C 2 that contain the static quarks and join the ends of the light quark 
propagator, as indicated in Fig. [|. 


I 11 eq (|28|) we have suppressed the spin labels of the quarks. However a feature of 
diquark dynamics is its sensitivity to the spins of the constituent quarks |19|] . For this 
reason it is important in this case to take account of the static quark spins as well as 
that of the light quark. We will do this by assuming that the static quark and light 
quark bind in a singlet state. The triplet state, being higher in mass |L9| , we will ignore 
for simplicity. In that case lines in our diagramatic model that carry a diquark state 
will have simple scalar propagators while the lines carrying a light or static quark in 
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the direction e will have the standard spin factor (1 + 7 .e )/2 as well as the appropriate 
energy exponential. 

A typical diagram that contributes to the amplitude in our simplified planar diagram 
model is shown in Fig. [|. The light quark propagates along the static quark lines 
in a diquark combination and also makes transitions between the static quarks on 
each side of the diagram. This diagram is part of the quenched approximation. The 
incorporation of dynamical quarks into the model results in diagrams such as that in 
Fig. 1 where internal quark loops appear creating intermediate states that involve a 
static nucleon, containing one static and two light quarks, together with a static meson, 
containing a static quark and a light anti-quark. 

7.1 Quenched Calculation 

Because the light quark can make a transtion from one static quark to the other, 
even the quenched calculation is not trivial. Within the model we are concerned with 
two channels. The left channel in which the light quark propagates along the static 
quark at the origin and the right channel in which it propagates along the static quark 
displaced a distance R from the origin. We will associate a hopping parameter k" with 
the motion of the light quark along a static quark line. We will introduce an energy 
Wd that provides a factor yJWd for the strength of emission of the light quark from a 
static diquark. With these rules the contribution to J-{T) of the diagram in Fig. [4] is 

(2 K"e~ rTR ) Tl (^V m ^) (2K"e~ aR ) T2 ( 2K"e ~ aR ) J3 (29) 

2 y 2 V 2 y 

The factors of ~ that accompany the passage of the light quark across the diagram 
follow from the 7 -matrix algebra. Because it plays no further role from now on we will 
drop the quark spin factor in T{T) . In order to sum up diagrams of the type we have 
just evaluated we require a 2 x 2 matrix of propagators in order to describe the two 
channels in the problem. We have 

Qt T \ _ ( g ll(T) G lr (T) \ 

m ~ V G ^(T) G rr (T) ) ■ (30) 

Up tp terms that vanish with the lattice spacing, J-(T) = Gll{T) . The summation 
over diagrams proceeds by requiring the propagators to obey the stepping relation 

G(T + 1) = A'G(T) , (31) 
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where 


(32) 


A! 


a' a'd \ 

a'd a' J ’ 


where a' = 2n"e~ aR = e^ cR ~ Ed and c = \Wde~ mqR and we have identified Ed = 
— log(2/c") as the static diquark energy. It is immediately obvious that there are two 
eigenchannels. A symmetrical superposition of the left and right channels and an 
antisymmetric superposition. 

Choosing these as the basis channels A' takes the form 


A' 


a'( 1 + d) 

a'( 1 



(33) 


The eigenenergies of these channels respectively e± = — log(a'(l ± d)) ~ Ed + aR ± 
\Wde~ mqR . We see here the influence of the quark exchange in splitting the degeneracy 
of the two string-diquark channels. 


7.2 Dynamical Quarks 

The effect of dynamical quarks is to insert light quark loops into the diagrams. As 
pointed out above this creates the possibility of two new channels one with a static 
nucleon at the origin and a static meson at a distance R and the other with the nucleon 
and meson interchanged. See Fig [5] . 

In order to sum over the diagrams of the model we need a 4 x 4 matrix of correlators. 


G(T) = (Gy(T)} 


(34) 


where i,j = 1,2, 3,4 and the labels 1 and 2 refer to the diquark-string channels with 
the diquark on the left and right respectively (L and R above). The labels 3 and 4 
refer to the nucleon-meson channels with the nucleon on the left and right respectively. 
The upgrade step is 


where now 


G(T + 1) = PG(T) , 


/ 

a' 

a'd 

a'd' 

0 

\ 


a'd 

a' 

0 

a'd' 


b'd' 

0 

b' 

0 


V 

0 

b'd' 

0 

b' 

J 


where a' and d are as before and 

b' = e - (E » +E » ) , and c" 


VWWR- m ,n 

V6 


(35) 


(36) 


(37) 
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Here Wn is the energy parameter associated wiith the emission of a light quark by a 
static nucleon. The mixing phenomena in the model are easily disentangled because the 
left-right symmetry of P means that the symmetric diquark-string state mixes only with 
the symmetric nucleon-meson state and similarly for the corresponding antisymmetric 
states. If we choose these symmetric and anti-symmetric states as the basis channels 
then G{T) and P have the form 

G^ S \T) 0 \ p_(P {S) 0 \ 

o gM\t) ) ’ ^ “ v 0 p{A) ) ’ 




and 


p(S) _ ( + d) a'c" \ p (A) 

^ b'c" b' ) 1 


a!{l — d) a'c" 
b'c" b' 


(39) 


We have G (S) (T + 1) = P^G (S) (T) and G ( - a) (T+ 1) = P^G (A) {T) with the solutions 
G (S)( T ) = (p(S)^j T and q(A)^ = (p(A)y We can write 


p(S,A) — PQ{S,A)^(S,A)pq(S,A)^P~^ 


(40) 


where 


f Va 7 0 \ q(S,a) = f co s &s,a ~ sin 9 s .a \ 
\ 0 \[b' ) ’ \ sin 9 s ,a cos 6 s ,a ) 


(41) 


The eigenvalues of pGd) 
in the diagonal matrix 


are the columns of DO^ s - A ) and the eigenvalues are the entries 


\(S,A) 


xf A) 0 \ 

0 aL s ’ a) J 


(42) 


It is easily established that 


i |a'(l + d) + b' ± ^(a'(l+c') - b') 2 + Aa'b'c" 2 ^ (43) 

and 

\ { ± ] = i ja'(l - d) + b'± y/(a'( 1 - d) - 6') 2 + 4a/6'c //2 | (44) 

The mixing angles are given by 


tan 6 s 


—(a'(l + d) - b') + (a'(l + d) - b') 2 + Aa'b'c" 2 

2 yfdVd' 


(45) 


and a similar formula for tan 6a where d —> —d . The mixing analysis is very similar to 
that for the original string breaking model and leads to the results that the critical value 
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of the separation in the symmetric and antisymmetric cases are given by b' — a(l ± c'), 
that is 

oR c = E n + E m - E d ± Uv der rn q R c ( 46 ) 

and the mixing ranges are 


A R s ,a = vr 


[2 yww^e~ m ^ 

V 3 cr ± \m q Wde~ rn i Rc 


and the energy splits at maximum mixing are 


(47) 


A E s ,a = ^\\lwW~Ne- m " R ‘ . (48) 

Although the energies are likelier to be higher than in the standard string scenario 
the static quark-diquark string is also of great interest as an alternative string system 
specific to SU (3) and because of what it can reveal about diquark dynamics. 


8 Conclusions 

We have explored further the consequences of a simple picture of string breaking as a 
mixing phenomenon. We have clarified the elimination of lattice spacing dependence 
in the model by improving it over its original formulation with the introduction of 
appropriate energy factors for quark and meson emission from bound states. We have 
shown that it is possible to incorporate a meson-meson potential into the the original 
model and find, as might be expected, that it does not lead to radically different results 
on mixing. 

We have also investigated the effect of degenerate light quark flavours. Here the only 
significant effect is a slight strengthening of the mixing effect by a factor \JN/ . It 
is also interesting that only the totally symmetric combination of two-meson channels 
participates in string breaking leaving the orthogonal combinations unmixed. 

We point out that a comparison of the quenched and unquenched calculations of string 
breaking is of considerable interest as a simulation. The appropriate measurement is 
of the string-meson transition amplitude. On the basis of the model we show that the 
two cases are likely to be rather similar in the outer part of the mixing region but 
show distinct behaviours in the central part of the mixing region, where the quenched 
amplitude shows a sharp transition between energy exponentials and the unquenched 
case shows a more gradual crossover. It is possible, if the comparison between the 
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two cases is made in an appropriate way in an actual simulation, that the quenched 
measurement could provide a prediction for the energy split and the size of the mixing 
region in the unquenched case. 

We study static quark-diquark string breaking which shows some simple but interesting 
patterns of channel mixing. The basic phenomenon is very similar to that of the quark- 
anti-quark string breaking. The symmetry of the system however provides us with two 
channels in the quenched case. These each go on to mix with appropriately symmetrical 
or antisymmetrical meson-nucleon channels. This pattern of double string breaking is 
interesting in itself and more so since it is specific to the case of SU(3) . It seems worth 
investigating this case even though the higher energies involved will make it harder to 
measure. Perhaps an approach using a non-symmetric lattice with a fine lattice spacing 
the time direction would be appropriate in this case. 

The success of simulations using the SU(2) Higgs model in observing string breaking 
H, U suggests that similar results could be obtained for the scenarios discussed in this 
paper with scalar matter fields replacing the dynamical quarks. It would of course be 
necessary to go to the SU( 3) Higgs model to see the effects discussed in the previous 
section. 
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Figure 1: Wilson loop (h 


taining internal quark loop (light line). 
























